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Abstract. This paper is concerned with the recognition of approximate graph products with re- 
spect to the Cartesian product. Most graphs are prime, although they can have a rich product-like 
structure. The proposed algorithms are based on a local approach that covers a graph by small sub- 
graphs, so-called partial star products, and then utilizes this information to derive the global factors 
and an embedding of the graph under investigation into Cartesian product graphs. 
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1. Introduction 

This contribution is concerned with the recognition of approximate products with respect to the 
Cartesian product. It is well-known that graphs with a non-trivial product structure can be recog- 
nized in linear time in the number of edges for Cartesian product graphs [16]. Unfortunately, the 
application of the "classical" factorization algorithms is strictly limited, since almost all graphs are 
prime, i.e., they do not have a non-trivial product structure although they can have a product-like 
structure. In fact, even a very small perturbation, such as the deletion or insertion of a single edge, 
can destroy the product structure completely, modifying a product graph to a prime graph [3, 22]. 
Hence, an often appearing problem can be formulated as follows: For a given graph G that has a 
product-like structure, the task is to find a graph H that is a non-trivial product and a good approxi- 
mation of G, in the sense that H can be reached from G by a small number of additions or deletions 
of edges and vertices. The graph G is also called approximate product graph. 

The recognition of approximate products has been investigated by several authors, see e.g. 
[4, 9, 10, 8, 17, 22, 14, 20, 21, 6, 11]. In [17] and [22] the authors showed that Cartesian and strong 
product graphs can be uniquely reconstructed from each of its one-vertex-deleted subgraphs. More- 
over, in [ 18] it is shown that fc-vertex-deleted Cartesian product graphs can be uniquely reconstructed 
if they have at least k+ 1 factors and each factor has more than k vertices. In [14, 20, 21] algorithms 
for the recognition of so-called graph bundles are provided. Graph bundles generalize the notion of 
graph products and can also be considered as a special class of approximate products. Equivalence 
relations on the edge set of a graph G that satisfy restrictive conditions on chordless squares play a 
crucial role in the theory of Cartesian graph products and graph bundles. In [11] the authors showed 
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that such relations in a natural way induce equitable partitions on the vertex set of G, which in turn 
give rise to quotient graphs that can have a rich product structure even if G itself is prime. However, 
Feigenbaum and Haddad proved that the following problem is NP-complete 

Problem 1.1 ( [4]). To a given connected prime graph G find a connected Cartesian product 
G\ □ . . . OGk with the same number of vertices as G, such that G can be obtained from G\ □ . . . QG^ 
by adding a minimum number of edges only or deleting a minimum number of edges only. 

Hence, in order to solve this problem not only for special classes of graphs but also for gen- 
eral cases one should provide heuristics that can be used in order to solve the problem of finding 
"optimal" approximate products. A systematic investigation into approximate product graphs w.r.t. 
the strong product showed that a practically viable approach can be based on local factorization 
algorithms, that cover a graph by factorizable small patches and attempt to stepwisely extend re- 
gions with product structures [9, 10, 8]. In the case of strong product graphs, one benefits from the 
fact that the local product structure of induced neighborhoods is a refinement of the global factors 
[8]. However, the problem of finding factorizable small patches in Cartesian products becomes a bit 
more complicated, since induced neighborhoods are not factorizable in general. In order to develop a 
heuristic, based on factorizable subgraphs and local coverings which in turn can be used to factorize 
large parts of the possibly disturbed graph we introduce the so-called partial star product (PSP). The 
partial star product is, besides trivial cases such as squares, one of the smallest non-trivial subgraphs 
that can be isometrically embedded into the product of so-called stars, even if the respective induced 
neighborhoods are prime. Considering a subset of the set of all partial star products of a graph, we 
propose in this contribution several algorithms to compute so-called product colorings and coordi- 
natizations of the subgraph induced by the partial star products. This information can then be used 
to embed large factorizable subgraphs of possibly prime graphs into a Cartesian product. 

This contribution is organized as follows. We begin with an introduction into necessary pre- 
liminaries and continue to define the partial star product. We proceed to give basic properties of the 
partial star product and concepts of product relations based on PSP's. These results are then used to 
develop algorithms and heuristics that compute (partial) factorizations of given (un)disturbed graphs. 



2. Preliminaries 
2.1. Basic Notation 

We consider finite, simple, connected and undirected graphs G=(V,E) with vertex set V (G) = V and 
edge set E(G) = E. A map y:V(H) ->■ V(G) such that (x,y) G E(H) implies (y(x),y(y)) e E(G) 
for all x,y £ V(G) is a homomorphism or embedding of H into G. We call two graphs G and H 
isomorphic, and write G ~ H, if there exists a bijective homomorphism y whose inverse function is 
also a homomorphism. Such a map y is called an isomorphism. 

For two graphs G and H we write GUH for the graph (V (G)UV (H),E(G)0E(H)), where U 
denotes the disjoint union. The distance dc(x,y) in G is defined as the number of edges of a shortest 
path connecting the two vertices x,y <G V(G). A graph H is a subgraph of a graph G, in symbols 
H C G, if V(H) C V(G) andE(tf) C E(G). A subgraph// C G is isometric i£d H (x,y) = d G {x,y) for 
all x,y € V(H). For given graphs G and H the mapping y : V(H) — > V(G) is an isometric embedding 
if d(j(u,v) = da (/(«), 7(v)) for all u,v € V(G). For simplicity, in such case we also call H isometric 
subgraph of G. If H C G and all pairs of adjacent vertices in G are also adjacent in H then H is called 
an induced subgraph. The subgraph of a graph G that is induced by a vertex set W C V (G) is denoted 
by (W). An induced cycle on four vertices is called chordless square. Let the edges e = (v,u) and 
/ = (v, w) span a chordless square ({v,u,x, w}). Then / is the opposite edge of (x,u). The vertex x is 
called top vertex (w.r.t. the square spanned by e and /). A top vertex x is unique if \N[x] PliV[v] | = 2. 
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In other words, a top vertex x is not unique if there are further squares with top vertex x spanned by 
the edges e or / together with a third distinct edge g. 

We define the open k-neighborhood of a vertex v as the set A^(v) = {x G V(G) | < dc(v,x) < 
k). The closed k-neighborhood is defined as A^[v] = A^(v) U {v}. Unless there is a risk of confusion, 
an open or closed k-neighborhood is just called k-neighborhood and a 1 -neighborhood just neighbor- 
hood and we write N(v), resp. N[v] instead of Ni (v), resp. N\ [v]. To avoid ambiguity, we sometimes 
write N^(v), resp. N^[v] to indicate that Nk(v), resp. A^[v] is taken with respect to G. 

The degree of a vertex v is defined as the cardinality |iV(v)|. A star G = (V,E) is a connected 
acyclic graph such that there is a vertex x that has degree |V| — 1 and the other |V| — 1 vertices have 
degree 1 . We call x the star-center of G. 

2.2. Product and Approximate Product Graphs 

The Cartesian product GUH has vertex set V(GUH) = V{G) x V(H); two vertices (gi,h), {gi,h 2 ) 
are adjacent in GUH if {g\,g-i) € E(G) and h\ = ho, or (Ai,/^) € E{G-i) and gi = g2- The one-vertex 
complete graph /T] serves as a unit, as K\ UH = H for all graphs H. A Cartesian product GUH is 
called trivial if G ~ or // ~ Ai. A graph G is prime with respect to the Cartesian product if 
it has only a trivial Cartesian product representation. A representation of a graph G as a product 
G\ □ G% □ • • • □ Gfc of prime graphs is called a prime factor decomposition (PFD) of G. 

Theorem 2.1 ([19, 16]). Any finite connected graph G has a unique PFD with respect to the Carte- 
sian product up to the order and isomorphisms of the factors. The PFD can be computed in linear 
time in the number of edges of G. 

The Cartesian product is commutative and associative. It is well-known that a vertex x of a 
Cartesian product U" =1 Gj is properly "coordinatized" by the vector c(x) := (c\ (x), . . . ,c„(x)) whose 
entries are the vertices c, (x) of its factor graphs G, [7]. Two adjacent vertices in a Cartesian product 
graph therefore differ in exactly one coordinate. Note, the coordinatization of a product is equivalent 
to an edge coloring of G in which edges (x,y) share the same color c^ if x and y differ in the 
coordinate k. This colors the edges of G (with respect to the given product representation). It follows 
that for each color c the set E c — {e £E(G) \ c(e) = c} of edges with color c spans G. The connected 
components of (E c ), usually called the layers or fibers of G, are isomorphic subgraphs of G. A partial 
product H C G is an isometric subgraph of a (not necessarily non-trivial) Cartesian product graph 
G. 

For later reference, we state the next two well-known lemmas. 

Lemma 2.2 (Distance Lemma, [12]). Let x = (xq,xh) andy = (yG,yti) be arbitrary vertices of the 
Cartesian product of GUH. Then 

dcoH{x,y) = d G {x G ,y G ) + d H {x H ,y H ) . 

Lemma 2.3 (Square Property, [12]). Let G = D" =1 G/ be a Cartesian product graph ande,f GE(G) 
be two incident edges that are in different fibers. Then there is exactly one chordless square in G 
containing both e and f. 

For more detailed information about product graphs we refer the interested reader also to [7, 
12] or [13]. 

For the definition of approximate graph products we defined in [9] the distance d(G,H) be- 
tween two graphs G and H as the smallest integer k such that G and H have representations G', 
H' for which the sum of the symmetric differences between the vertex sets of the two graphs and 
between their edge sets is at most k. That is, if 

\V(G')AV(H')\ + \E{G')AE(H')\<k. 
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A graph G is a k-approximate graph product if there is a product H such that 

d(G,H)<k. 

Here k need not be constant, it can be a slowly growing function of |2?(G)|. Moreover, the next 
results illustrate the complexity of recognizing approximate graph products. 

Lemma 2.4 ([9]). For fixed k all Cartesian k-approximate graph products can be recognized in 
polynomial time in n. 

Without the restriction on k the problem of finding a product of closest distance to a given 
graph G is NP-complete for the Cartesian product [4]; see Problem 1.1. 

2.3. Relations 

We will consider equivalence relation R on edge sets E, i.e., R C E x E such that (i) (e,e) G R 
(reflexivity), (ii) (e,f) G R implies (f,e) G R (symmetry) and (iii) (e,f) G R and (f,g) G R implies 
(e,g) G R (transitivity). We will furthermore write <p c R to indicate that <p is an equivalence class of 
R. A relation Q is finer than a relation i? while the relation R is coarser than 2 if (e,/) £ G implies 
(e,f) G ii, i.e, Q CR.ln case, a given reflexive and symmetric relation R need not be transitive, we 
denote with R* its transitive closure, that is the finest equivalence relation on E(G) that contains R. 
For a given graph G = (V,E) and an equivalence relation R on E we define the R-coloring of G as a 
map of the edges onto its equivalence class, i.e, the edge e G E is assigned color k iff e G % E ii. 

For a given equivalence class <p c ^ and a vertex m G V(G) we denote the set of neighbors of u 
that are incident to u via an edge in <p by Nm(u), i.e., 

iV ? (»):={v£y(G) | [u,v]ep}. 

The closed (^-neighborhood is then Ny [u] = Ny (u) U {u}. 

For later reference we need the following simple lemma. 

Lemma 2.5. Let R be an equivalence relation defined on the edge set of a given graph G = {V,E) 
and H C G be a subgraph of G. Then the restriction R\ H = {(e,f) G R | e,f G E(H)} of R on the 
edge set E(H) is an equivalence relation. 

Proof. Clear. □ 

For the recognition of Cartesian products the relation 8 is of particular interest. 

Definition 2.6. Two edges e,f G E(G) are in the relation 8(G), if one of the following conditions 
in G is satisfied: 

(i) e and / are adjacent and there is no unique chordless square spanned by e and /. 

(ii) e and / are opposite edges of a chordless square. 

(iii) e=f. 

If there is no risk of confusion we write 8 instead of 8(G). Clearly, the relation 8 is reflexive 
and symmetric but not necessarily transitive. However, the transitive closure 8* is an equivalence 
relation on E(G) that contains 8. Note, that our definition of 8 slightly differs from the usual one, 
see e.g. [18, 15], which is defined analogously without forcing the cordless square in Condition 
(i) to be unique. However, for our purposes this definition is more convenient and suitable to find 
the necessary local information that we use to define those factorizable small patches which are 
needed to cover the graphs under investigation and to compute the PFD or approximations of it with 
respect to the Cartesian product. Moreover, as stated in [18, 15], any pair of adjacent edges that 
belong to different 8* classes span a unique chordless square, where 8 is defined without claiming 
"uniqueness" in Condition (/). Thus, we can easily conclude that the transitive closure of our relation 
8 and the usual one are identical. 
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Finally, two edges e and / are in relation o(G) if they have the same Cartesian colors with 
respect to the prime factorization of G. We call <j(G) the product relation. The first polynomial 
time algorithm to compute the factorization of a graph explicitly constructs a starting from the finer 
relation 8 [5]. The product relation a was later shown to be simply the convex hull €(8) of the 
relation 8(G) [15]. Notice that 8(G) C 8(G)* C a(G) [15]. 

3. The Partial Star Product 
3.1. Basics 

In order to compute 8 from local coverings of the graph G = (V,E) we need some new notions. 
Clearly, 8 is still defined in a local manner since only the (non-)existence of squares are consid- 
ered and thus, only the induced 2-neighborhoods are of central role. However, although the 2- 
neighborhood can be prime, we define subgraphs of 2-neighborhoods, that are factorizable or at 
least graphs that can be isometrically embedded into Cartesian products and have therefore a rich 
product structure. For this purpose we define for a vertex v G V(G) the relation d v , that is a subset 
of 8 and provides the desired information of the local product structure of the subgraph (A^[v]). 
Based on the transitive closure 0* we then define the so-called partial star product S v , a subgraph 
of (N2 [v]), which provides the details which parts of the induced 2-neighborhood are factorizable or 
can be isometrically embedded into a Cartesian product. 

Let G — (V,E) be a given graph, v G V and E v be the set of edges incident to v. The local 
relation d v is then defined as 

V = ((E v xE)U(Ex E v )) n 8(G) C 5 ( (JV 2 G [v] } ) . 

In other words, V is the subset of 8(G) that contains all pairs (e,f) G 8(G), where at least one of 
the edges e and / is incident to v. Note, 0* is not necessarily a subset of 8 but it is contained in 8*. 
For a subset W C V we write Q V (W) for the union of local relations V , v G W: 

d v (W) = U veW X) v . 

We now define the so-called partial star product S v , that is, a subgraph containing all edges 
incident to v and all squares spanned by edges e,e' G E v where e and e' are not in relation 0*. To be 
more precise: 

Definition 3.1 (Partial Star Product (PSP)). Let F r CE\E V be the set of edges which are opposite 
edges of (chordless) squares spanned by e,e' G E v that are in different D* classes, i.e., (e,e') £ 0*. 

The partial star product is the subgraph S v C G with edge set E 1 = E v U F v and vertex set 
\J eeE ie. We call v the center of S r , edges in E v primal edges, edges in F r non-primal edges, and the 
vertices adjacent to v primal vertices with respect to S v . 

The reason why we call S v a partial star product is that S v is an isometric subgraph or even 
isomorphic to a Cartesian product graph H of stars, as we shall see later (Theorem 3.9). Hence, S v 
is a partial product of H. For the construction of this graph H we introduce the so-called star factors 
Si, see also Figures 1 and 3. 

Definition 3.2 (Star Factor). Let G = (V,E) be an arbitrary given graph and S v be a PSP for some 
vertex v G V. Assume 9* has equivalence classes tfa, . . . , (p n . We define the star factor S; as the graph 
with vertex set [v] that contains all primal edges of E v that are also in the induced closed <p,- 
neighborhood, i.e., £(§,•) = E^N^ [v]))n£ v . 

Note, this definition forbids triangles in and hence, each S, is indeed a star. We denote the 
restriction of 5* to the subgraph 5 V with 

*\ Sr :=K\ Sv = {(eJ)£K\ej£E(S v )}. 
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Figure 1 . Examples of various PSP's S v highlighted by thick edges. Note, in all 
cases except in case (/) the set F v is empty and hence, the PSP's S v in the other 
cases just contain the edges incident to v. 




Figure 2. Left: A hypercube Q3 is shown. The three equivalence classes of 
8* (Qi) are highlighted by solid, dashed and double lined edges, respectively. 
Right: The PSP S v is shown. Again, has three equivalence classes. How- 
ever, since the edges (0,1) and (1,2) as well as the edges (2,3) and (3,4) span 
no square we can conclude that 8* (S v ) just contains one equivalence class. Hence, 

In other words, 5i Sv is the subset of 0* that contains all pairs of edges (<?,/) € 3* where both edges 
e and / are contained in S v . We want to emphasize that V*, s ^ 8*(S V ); see Figure 2. In addition, by 
Lemma 2.5 we can conclude that 0i Sv is an equivalence relation. For a given subset W C V we define 

5| Sv (W) = U v€W t)| Sv 

as the union of relations 0| Sv , v £ W. As it will turn out, for a given graph G = (V,E) the transitive 
closure 5|s,,(^)* is the equivalence relation 8(G)*, see Theorem 3.11. 

3.2. Properties of the Partial Star Product 

We now establish basic properties of the graph 5V, its edge sets E v and F v , as well as of the relation 
5* and its restriction 0i Sv to S v . 

Lemma 3.3. Given a graph G = (V,E) and a vertex v £ V. Then F v = if and only if for all edges 
e,e' £ E v holds (e,e') £ 0*. Moreover, ifF v ^ then \F V \ > 2. 

Proof. Clearly, if for all edges e, e' £ E v holds (e, e') £ 0* then by definition F v = 0. 
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Figure 3 . Shown is a graph G ~ (N? [v] ) . Note, 8(G)* has one equivalence class 
and thus, G is prime. However, the partial star product (PSP) S v , that is the sub- 
graph that consists of thick and dashed edges is not prime. The subgraph S v is 
isomorphic to the Cartesian Product of a star with four and a star with three ver- 
tices. The two equivalence classes of V\g v are highlighted by thick, resp. dashed 
edges. 



Let F v — and assume there are edges e,e' <E E v that are not in relation 0*. In particular, these 
edges are not in relation 0,,, and therefore not in relation 8(G). By Condition (/) of Def. 2.6 and since 
e and e' are adjacent, there is a chordless square containing e and e' and therefore, respective opposite 
edges / and /'. Condition (ii) of Def. 2.6 implies (e,f),(e',f') G 8(G). Therefore, /,/' G F v , a 
contradiction. 

Furthermore, since F v contains all opposite edges of squares spanned by e,e' G E v we can easily 
conclude that \F V \ > 2, if F v ^ 0. □ 

Lemma 3.4. Let G=(V,E) be a given graph and let S v be a PSP for some vertex v G V. If e,f G E v 
are primal edges that are not in relation t)*„ then e and f span a unique chordless square with a 
unique top vertex in G. 

Conversely, suppose that x is a non-primal vertex ofS v , then there is a unique chordless square 
in S v that contains vertex x and that is spanned by edges e,/ G E v with (e,f) G" 5*. 

Proof. First, we show that e and / span a unique chordless square. By contraposition, assume e and 
/ span no unique chordless square. Since e and / are adjacent, Condition (;) of Def. 2.6 implies 
that (e,f) G 8(G) and hence, (e,f) G 5 V Q D*. Therefore, if (e,f) £ 0*, then they must span a 
unique chordless square. Let e = (v, u) and / = (v, w), (e,f) ^ t>*, span the unique chordless square 
SQi = ({v,",!^}) and assume for contradiction that the top vertex x is not unique. Hence, there 
must be at least three squares: the square SQi, the square SQ2 = ({v, u,x,y}) spanned by e and 
g, and the square 5^3 = ({v, w,x,y}) spanned by / and g = (v,y). We denote edges as follows: 
a = (x,y) and b = (x, w). Assume both squares SQ2 and SQj are chordless. Then Def. 2.6 (ii) implies 
(f,a),(a,e) G 8(G) and therefore, (e,f) G 0*, a contradiction. If both squares have a chord then 
Def. 2.6 (i) implies that (e,g),(f,g) G 8(G) and thus, (e,f) G U*, again a contradiction. If only one 
square, say SQ2, has a chord (u,y), then (e,g) G 8(G) and (f,a), (g,a) G 8(G) and again we have 
(e,f) G K- 

Assume x is a non-primal vertex in S v . By definition, there are non-primal edges /' = (x, u) , e' = 
(x,w) G F v that are contained in a square spanned by e — (v,u),f = (v,w) G E v , whereas (e,/) ^ f*. 
As shown above, the square spanned by e and / is unique with unique top vertex. Hence, if there is 
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another square in S v containing x then it must be spanned by e' , f and this square contains additional 
edges /" = (y, u) , e" — (y, w). However, then there is a square ({v, u,y, w}), which contradicts the fact 
that the square spanned by e and / is unique. If the unique square spanned by e and / is not chordless, 
then Def. 2.6 (?) implies (e,f) G S(G) and thus (e,f) G Z>*, a contradiction. □ 

By means of Lemma 3.3 and 3.4 and the definition of partial star products we can directly infer 
the next corollary. 

Corollary 3.5. Let G=(V,E) be a given graph and let S v be a PSP for some vertex v G V. 

1. If (e,f) G D* then there is no square in S v spanned by e and f. 

2. Every square in S v contains two edges e, e' G E v and two edges /, /' G F v , and every edge f G F v 
is opposite to some primal edge e G E v . 

3. Every non-primal vertex in S v is a unique top vertex of some square spanned by edges e,e' G E v . 

Lemma 3.6. Let G=(V,E) be a given graph and let f G F v be a non-primal edge of a PSP S v for 
some vertex v G V. Then f is opposite to exactly one primal edge e G E v in S r and (e,f) G d\g v . 

Proof. By Corollary 3.5, construction of S v and since / G F v , there is at least one edge e G E v such 
that / is opposite to e and therefore at least one square SQ\ = ({v,w,x, u}) in S v spanned by primal 
edges e = (v,«) and e' = (v,w) that contains the edge / = (w,x). Note, by construction (e,e') g" D* 
and e is opposite to /. Assume for contradiction that / is opposite to another edge g = (v,y). Then 
there is another SQ2 = ({v,y,x, w}). Hence, e and e' do not span a square with unique top vertex in 
G. By Definition 2.6 and Lemma 3.4 we can conclude that (e,e') G 3*, a contradiction. Hence e and 
e' span a unique chordless square containing the edge /. By Condition (i) of Definition 2.6 it holds 
(e,f) G 8. Since e G E v we claim (<?,/) G d v and consequently (e,f) G f| s • □ 

Lemma 3.7. Let G=(V,E) be a given graph with bounded maximum degree A and ffCV such that 
(W) is connected. Then each vertex meets every equivalence class ofD^ v (W)* in U v ew^v. i- e -> 

for each equivalence class (p E f|s v (W)* and for each vertex j£ff there is an edge (x,y) G <p with 
(x,y) G E(U v£ wS v ). Moreover, f (W)* has at most A equivalence classes. 

Proof. Let v G W be an arbitrary vertex and S v be its PSP. We show first that v meets every equiv- 
alence class of D|5 v in S v . Assume for contradiction that there is an equivalence class <p e 5^ that 
is not met by v and hence for all edges e G E v we have e £ (p. Hence, there must be a non-primal 
/ G F v with / G (p. By construction of S v and by Lemma 3.6 this edge / is opposite to exactly one 
edge e G E v with (e,f) G 0|s v , but then e G <p, a contradiction. We show now that every primal ver- 
tex w in S v meets every equivalence class of 0i$ v . Let (p e 0k v be an arbitrary equivalence class. If 
e = (v,w) G 9 we are done. Therefore assume e g" (p. Hence, there must be at least a second equiv- 
alence class <p' e 0i Sv with e G <f>'. Since vertex v meets every equivalence class there is an edge 
e' = (v,u) G (p. Moreover, since (e,e') 0* it follows that (e,e') V C 5. Since e and e' are adja- 
cent and by Condition (i) of Definition 2.6 the edges e and e' span a unique chordless square. Hence, 
there is an opposite edge / = (w,x) of e' . By construction of S v we have / G F v and hence, Lemma 
3.6 implies (e',f) G 0|5 V . Therefore, the primal vertex w meets equivalence class <p in S v . Note, not 
every equivalence class of must be met by non-primal vertices in S v in general, as one can easily 
verify by the example in Figure 4. 

It remains to show that every vertex x G W meets every equivalence class of 5|s v (W)* in 
UvewSv Assume we have chosen an arbitrary vertex x G W, computed S x and 0| St . As shown, vertex 
x and all its primal neighbors y in S x meet every equivalence class of t)is x - Assume W contains more 
then one vertex. Since (W) is connected there is a primal vertex y of x that is also contained in W. 
Hence, vertex x is a primal neighbor of y in S y and every equivalence class of Di s is met by y as well 
as by x. Let (p e (3^^ U5|s y )* be an arbitrary equivalence class. Assume neither x nor y meets (p. Then 
each edge / G <p must be in or /y Assume / G F y then, by construction of Sy and Lemma 3.6, this 
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edge / is opposite to exactly one edge e = (y,a) G E y with (e,f) G , and hence e G <p, a contra- 
diction. Assume now all edges e G 9 are only met by y but not by and therefore, e' = (x,y) g" 9. 
However, since e and e' are in different equivalence classes of (0| 5r U0|j v )* they must be in dif- 
ferent equivalence classes of T)w . Hence, (e,e') ^ 0* and thus, (e,e') £ V C 8. Since e and <?' are 
adjacent and, by Condition (;) of Definition 2.6, the edges e and e' span a unique chordless square. 
Hence, there is an opposite edge / = (x,w) of e in S y and, by Lemma 3.6 we conclude (e,f) 6 O15 
and therefore, f € (p, which implies that x meets 9, a contradiction. Hence, every equivalence class 
(p E (0| Sr UDi 5y )* must be met by x and y. By the same arguments one shows that each primal vertex 
of x in S x and y in 5), meets every equivalence class of (0| St U 5|s v )*- If W \ {x, y} / we we can 
choose a primal neighbor z G W of x or y, since (W) is connected. By the same arguments as before, 
one shows that each vertex x,y, resp. z and each of its primal vertices in S x ,S y , resp. S z meets every 
equivalence class of ((D|s t U?)^ )* U0|s_)* = (V\s x U5|s v U0|5. )* in 5 X U S y U S,. Therefore, we can 
traverse (W) in breadth-first search order and inductively conclude that every vertex x G W meets 
every equivalence class of i)|s v (W)* in U ve wSv 

Finally, we observe that each edge in E v might define one equivalence class of 0| Sv for each 
vertex v G W. Thus, 0js v can have at most A equivalence classes. Since this holds for all vertices and 
since equivalence classes in are combined equivalence classes of the respective 3| Sv classes, 

the number of equivalence classes in can not exceed A. □ 

In order to prove that each PSP can be isometrically embedded into a Cartesian product of 
stars, which is shown in the next theorem, we first need the following lemma. 

Lemma 3.8. Let G = O l j=l Gj be the Cartesian product of the stars. Assume the vertices in each 
V(Gi) are labeled from 0, . . . , |V(G,)| — 1, where the vertex with label always denotes the star- 
center of each G;. Let vq be the vertex with coordinates c(vq) = (0, ... ,0) Then the induced closed 
k-neighborhood {N^[vg\} is an isometric subgraph of G. 

Proof. Let (N^ [vq] ) be the induced closed ^-neighborhood of vq in G. Let x, y G [vg] be arbitrary 
vertices and let / C { 1 ,...,/} be the set of positions where x and y differ in their coordinate. More- 
over, let Iq C / be the set of positions where either x or y has coordinate 0. By the Distance Lemma 
we have d G (x,y) = Z ie i 1 + Le/\/ 2. 

We now construct a path from x to y that is entirely contained in [vq] and show that this path 
is a shortest path. Let i G Iq and w.l.o.g. assume c, (x) = 0, otherwise we would interchange the role 
of x and y. By definition of the Cartesian product there is a vertex y' that is adjacent to vertex y with 
Cj(y') = Cj(y) for all j ^ i and c,(y') = 0. By the Distance Lemma, we have dcj(cj(vG),Cj(y)) — 
dGj{cj(v G ),Cj{y')) for all j ^ i and d Gi {ci(v G ),Ci(y)) = d Gi (0,a(y)) = 1 and d Gi (ci(v G ),Ci(y')) =0 
and thus, d G (v Gl y') < d G (v G ,y) < k, which implies thaty' G N^[v G ]. We add (y,y') to the pathf(jc,y) 
from x to y and repeat to construct parts of the path from x to y' in the same way until all i G Iq are 
processed. In this way, we constructed subpaths P(x,v) and P{w,y) of P(x,y), both of which are 
entirely contained in (N^[v G ]) and |f(x,v)| + |f(w,y)| = |/o|- We are left to construct a path from 
v to w that is entirely contained in NP[v]. Note that by construction v and w differ only in the i-th 
position of their coordinates where i G I\Iq and c,(v) — cj(x) = cj(y) = cj(w) for all j $ I\Iq, 
By the definition of the Cartesian product for each / G / \/o there are edges (v,v'), resp. (v',v") 
such that v,v' and v" differ only in the i-th position of their coordinates. Since ^ ~ c,(v) 
and by definition of the Cartesian product it follows that c, (v') = and v" can be chosen such that 
c;(v") = c,(y) = Cj(w) 7^ 0. By the Distance Lemma and the same arguments as used before it holds 
dG(v G ,v') = d G (v G ,v") — 1 = d G (v G ,v) — \<k and hence, v',v" G N^[v G ]. Therefore we add the 
edges (v, v'), resp. (v', v") to the path from x to y, remove i from 7\/o and repeat this construction 
for a path from v" to w until I \ Iq is empty. 
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FIGURE 4. Shown is a graph G ~ (A^H). Note, 5(G)* has one equivalence class. 
The partial star product (PSP) S v is the subgraph that consists of thick, double- 
lined and dashed edges. Moreover, S v can be isometrically embedded into the 
Cartesian Product of a star with two and two stars with three vertices. The three 
equivalence classes of are highlighted by thick, double-lined, resp. dashed 
edges. 

Hence we constructed a path of length \Iq\ + 2|/\/o| = Hiei 1 +L(G/\/ 2 = dc(x,y). Thus, this 
path is a shortest path from x to y. Since this construction can be done for any x,y G N^[vq] we can 
conclude that (N^[vg\) is an isometric subgraph of G. □ 

Theorem 3.9. Let G = (V,E) be an arbitrary given graph and S v be a PSP for some vertex v G V. 
Let H = rf =1 §; be the Cartesian product of the star factors as in Definition 3.2. Then it holds: 

(1) S v is an isometric subgraph of H and in particular, S v ~ (A^[(vi,...,Vjfc)]) where V; denotes 
the star-center of Si, i — 1 , . . . , k. 

(2) 0| S „ c 8(H)* c a{H). 

(3) The product relation o(H) has the same number of equivalence classes as 5|s y . 
Proof. Assertion (1): 

If 0* has only one equivalence class, then there is nothing to show, since S v ~ §i ~ H. Therefore, 
assume d*. has k > 2 equivalence classes. 

In the following we define a mapping y : V(S V ) —> V(H) and show that y is an isometric em- 
bedding. In particular we show that y is an isomorphism from S v to the 2-neighborhood (N^Ivh]) 
for a distinguished vertex vh G V(H). Lemma 3.8 implies then that this embedding is isometric. 

For a given equivalence class <p, e 0* let (v) = {vi , . . . , V;} be the <f> ( - -neighborhood of the 
center v and be the corresponding star factor with vertex set V (S,) = {0, 1 , . . . , /} and edges (Q,x) G 
E($j) for all (v, v x ) G S v . Let H = □*_!§, be the Cartesian product of the star factors. The center v of 
S v is mapped to the vertex Vh G V(H) with coordinates c(\>h) = (0, . . . ,0), the vertices vj G ^ (v) 
are mapped to the unique vertex u with coordinates c r (u) = for all r ^ i and c, (m) = j. Clearly, these 
vertices exist, due to the construction of Si,... and since V(H) = x* =1 V(Si). Note, that these 
vertices we mapped onto are entirely contained in the 1 -neighborhood N h [\'h] of vu- Now let x be a 
non-primal vertex in S v . Hence, by Lemma 3.4 and Corollary 3.5, there is a unique chordless square 
({v,Vi,x,Vj}) in S v with unique top vertex x. Thus, v,- and vj are the only common neighbors of x 
in S v . Moreover, by definition and Lemma 3.4, the edges (v,v,-) G (p r and (v,vj) G (p s are in different 
equivalence, i.e., c^s. Thus, we map x to the unique vertex u with coordinates c/(m) = for all 
I r,s and c r (u) = i and c s (u) = j. Again, this vertex exists, due to the construction of Si, ... 
and since V(H) = x* =1 V(S,-). This completes the construction of our mapping y. 
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We continue to show that the mapping y : V(S Y ) — > iV* [vh] is bijective. It is easy to see that by 
construction and the definition of the Cartesian product, each primal vertex x has a unique part- 
ner y(x) in iVj [vh] and vice versa. We show that this holds also for non-primal vertices in S v 
and vertices in N? [vh] \N^[vh]- First assume there are two non-primal vertices x and x' in S v 
that are mapped to the same vertex u in H. Thus, by construction of our mapping y, the vertex 
x' must have the same primal neighbors v,- and v, as x in S v . However, by Lemma 3.4 this con- 
tradicts that (v,v,-) G (p r and (v,Vj) G % span a unique square. Therefore, y is injective. Now, let 
u G N^Ivh] \ N^[vh] be an arbitrary vertex in H. By the Distance Lemma we can conclude that 
duivn^u) =L?=i^S,-(0j Ci{u)). Moreover, since dn(vH,u) =2andd$.(0,Cj(u)) < 1 forall/= l,...,k 
we can conclude that cIh(vh,u) = d§ r (0,c r (u)) +dg s (Q,c s (u)) for some distinct indices r and s. As- 
sume that c r (u) = i and c s (u) = j. By construction, the star factor § r contains the edge (0,;) and 
S. s the edge (0, j). Hence, there are edges e = (v,v,) £ (p r and / = (v,vy) G <p s in S v . Lemma 3.4 
implies that there is a unique chordless square spanned by e and / with unique top vertex y that is 
also contained in S v . By construction of y the vertex v is the unique vertex that is mapped to vertex u 
in H. Since this holds for all vertices u G [vh] \ [vh], and by the preceding arguments, we can 
conclude that the mapping y : S v — > [vh] we defined is bijective. 

It remains to show that y is an isomorphism from S v to [vh] ■ By construction, every primal 
edge (v,vj) G <p,- is mapped to the edge (v#,x), where x has coordinates c,(x) = for i ^ r and 
c r (x) = j. Hence, (v, Vj) G E v if and only if (y(v),y(v/)) G E((N2 [vh]))- Now suppose we have anon- 
primal edge (vj,y) G <p r . By Lemma 3.4, there is a unique chordless square with edges (v,vj) G <p,- 
and (v,vj) G <p s and hence, by construction of §, and § s and the definition of the Cartesian product, 
there are edges e = (vh,z) and / = (vh,z') in H where z differs from vh in the r-th position of its 
coordinate and z' differs from vjj in the s-th position of its coordinate. By the Square Property, there 
is unique chordless square in H spanned by e and / with top vertex y' that has coordinates c, (y' ) = 
for ; 7^ r,s, c r (y') = i ^ and c s (y') = j ^ 0- By the construction of y we see that (vy,y) G F v 
implies (y(vj),y(y)) = {z',y') G E ( (A^ [vh] ) ) ■ Using the same arguments, but starting from squares 
spanned by e = (\'h,z) and / = (v#,z') inH, one can easily derive that (z',/) G E((N2[\'h])) implies 

(7" 1 (z'),r-i(y')) = (vy,)')eF v . 

Finally, Lemma 3.8 implies that (Ni[\'h]} is an isometric subgraph of H and therefore, y : 
V(S V ) — > V(H) is an isometric embedding. 
Assertion (2) and (3): 

W.l.o.g. we treat the graph 5\, as a proper subgraph of H. Lemma 2.5 implies that 8 (S v )* C 8(H)*. We 
continue to show that d\ Sv = 0*,^ C 8(S V )*. Let v G V(G) be a center of the PSP SV, and // = D^LjS,-, 
where §,■ are the corresponding star factors (w.r.t. S v ). Let e,f G £(5',,) such that (<?,/) G 5|j r . By the 
definition of 0| Sv and S v , at least one of the edges e or / must be incident with v and (e,f) G5(G) C 
8(G)*. If e,f G £ v are both primal edges then Corollary 3.5 and (e,f) G £>|s v C D* imply that e and 
/ span no square in S v , and therefore (e,f) G 8(S V ) C 5(5',.)*. Suppose e G £ v is a primal edge and 
/ G F v is non-primal. By Lemma 3.6, / is opposite to exactly one primal edge e' where (/,e') G Z)|s v . 
If e = e', then e and / are opposite edges in a chordless square in S v and thus (e,f) G 5(5,.). If 
e ^ e', then (e,f), (f,e') G 0|s v implies that (e,e') G 0|s v and we can conclude from Corollary 3.5 
that there is no square spanned by e and <?' in S v . Hence, the square in S v containing e' and / must be 
spanned by e' and another edge e" G E v . Then this square SQ contains the edges e',e",f and some 
edge /' G F v . Moreover, since e' is the unique opposite edge of / and since (e,f) G D|s v C 8(G)*, 
the edges e and / must be adjacent, which is only possible if e = e" or e is a chord in SQ. If e — e", 
then this would contradict that there is no square spanned by e and e' . If e is a chord in SQ then this 
contradicts that e' and / are opposite edges in a chordless square in S v . Therefore e = e' and thus, 
(ej) G 8(S V ) C 8(S V )*. Consequently, t>\ Sv C 5(5,,)*. 
Hence, we have 

0|^=^| Sir C5(S v )*c5(H)*Cff(ff). 
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Moreover since the number of 0| Sr classes equals the number of prime factors of H it follows 
that and <j(H) have the same number of equivalence classes. □ 

By the construction of star factors, the Distance Lemma and Theorem 3.9, we can directly infer 
the next corollary. 

Corollary 3.10. Let G = (V,E) be an arbitrary given graph, S v be a PSPfor some vertex v G V and 
0* have k — 1 or 2 equivalence classes. Then 

We conclude this section with a last theorem which shows that the transitive closure of the 
union 0|s v (V) over all vertices and its relations d v , even restricted to S v , is 8(G)*. 

Theorem 3.11. Let G = (V,E) be a given graph andd\$ r (V) = U Ve vt> \s v - Then 

* lSv (vy = s(Gy. 

Proof. By definition d v C 8(G). Moreover, Lemma 2.5, resp. Theorem 3.9 imply that Di^ C 5* C 
8(G)* for all v G V(G). Thus, 0| Sv (7) C 8(G)*, and hence t ]Sv (V)* C 8(G)*. 

Let e,f G E(G) be edges that are in relation 8(G). By definition, (e,f) G D v for some v 6 V(G). 
If e = (u,v) and / = (w, v) are adjacent, then e and / are contained in the set E v of S v and therefore in 
f W v C 5(G)*. Assume, e = (m, v) and / = (x,y) are opposite edges of a chordless square containing 
the edges e,f and g = (v,x). For contradiction, assume (e,f) ^ 0|s v (V)* and hence (e,/) ^ |s v (V). 
Thus, for each v G V we have (e,f) Oi^ and therefore, by definition, there is no square spanned 
by edges e,e' G E v with (e,e') ^ 5* such that / is the opposite edge of e. In particular, this implies 
(e,g) G 5* and hence (e,g) G J)|s v . Analogously, one shows that G <)\s x . Since 0|5 v UD|5 r C 

0| Sv (V) we can infer that (e,f) G f|s v (V)*, a contradiction. □ 

Theorem 3.11 allows us to provide covering algorithms for the recognition of 8(G)* or of 
8(H)* for subgraphs H C G that are based only on coverings by partial star products. Note, if 
<j(G) = 8(G)*, then the covering of G by partial star products would also lead to a valid prime 
factorization. However, as most graphs are prime we will in the next section provide algorithms, 
based on factorizable parts, i.e., of coverings where the PSP's have more than one equivalence class 
5| 5r , which can be used to recognize approximate products. 

4. Recognition of Relations, Colorings and Embeddings into Cartesian 
Products 

In order to compute local colorings based on partial star products and to compute coordinates that 
respect this coloring we begin with algorithms for the recognition of 5 \s v (W)* and 8(G)*. 

Lemma 4.1. Given a graph G = (V,E) with bounded maximum degree A and a subset ffcy such 
that (W) is connected, then Algorithm 1 computes D|s v (W)* and U v ew^v in 0(|V|A 4 ) time. 

Proof. The Algorithm scans the vertices in an arbitrary order and computes (iVj [v]), 8' = 
5((A^t'[v]}), as well as S v and D| 5v w.r.t. 5'. In order to compute the transitive closure of t)\s v (W) 
an auxiliary graph, the color graph F, is introduced. For each vertex v and to each equivalence class 
of 0|5 v some unique color is assigned, and T keeps track of the "colors" of the equivalence classes. 
All vertices of T are pairs (e,c). Two vertices (e',c ! ) and (e",c") are connected by an edge if and 
only if there is an edge e G <jV H <p c " with <p c i e and (p c ii e for some u,w G W. In other 
words, if there is an edge e that obtained both, color c' and c". Edges in T "connect" edges of local 
equivalence classes that belong to the same global equivalence classes in f |s v (W)*. The connected 
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Algorithm 1 Local 0| Sv (W)* computation 

l: INPUT: A graph G = (V,E), W C V. 

2: a^W 

3: initialize graph T = 0; {called "color graph"} 

4: while a ^ do 

5: take any vertex v of <r; 

6: compute (JV 2 G [v]), 8' = 8({N?[v})), S v and J)* Sv w.r.t. <5'; 

7: color the edges of S v w.r.t. the equivalence classes of D| 5v ; 

8: set numjolass = the number of equivalence classes of 

9: add numjclass new vertices to T; 

10: for every edge e in S v do 

11: if e was already colored in G then 

12: x = old color of e\ y = new color of e; 

13: add vertices (x,e) and (y,e) to F 

14: join all vertices of the from (x,f) and (y,f') in T; 

15: end if 

16: end for 

17: delete v from a; 

18: end while 

19: {compute the equivalence class (jfy e D|s v (W)*.} 

20: set numjoomp = number of connected components of T; 

21: for k — 1 to numjcomp do 

22: if color of e is vertex in component k of T then 

23: (p k <— e; 

24: end if 

25: end for 

26: OUTPUT: 0| Sl ,0*T and U veW S v ; 



components Q of T define edge sets Eq = ^uaeqVc- We therefore can identify the transitive closure 
of fl|s v (W)* by defining e G (fig c D|5 v (W)* if e G £g. Finally, we observe that this is iteratively done 
for all vertices v GW, that all edges in E((W)) are contained in some E v of 5 V and, by Lemma 3.7, 
that every equivalence class of Q\s v (W)* is met by every vertex v GW. Therefore, we can conclude 
that each edge is uniquely assigned to some class <pg e d\ Sii (W)*. Hence, the algorithm is correct. 

In order to determine the time complexity we first consider line 6. The induced 2-neighborhood 
can be computed in A 2 time and has at most A 2 vertices, and hence at most A 4 edges. As shown by 
Chiba and Nishizeki [2] all triangles and all squares in a given graph G = (V,E) can be computed 
in 0(|£|A) time. Combining these results, we can conclude that all chordless squares can be listed 
in 0(\E\A) time. Thus, in this preprocessing step, we are able to determine 8' ,S V and d\g v in 0(A 4 ) 
time. Since this is done for all vertices v G W, we end in an overall time complexity 0(|£| A+ |W|A 4 ) 
for the preprocessing step and the while-loop. For the second part, we observe that F has at most 
0(|£|) connected components. Since the number of edges is bounded by | V | A we conclude that 
Algorithm 1 has time complexity 0(\V |A 2 + |W|A 4 ) = 0(|V|A 4 ). □ 

By means of Theorem 3.11 and Lemma 4. 1 we can directly infer the next corollary. 

Corollary 4.2. Let G = ( V, E ) be a given graph with bounded maximum degree A. Then 8(G)* can 
be computed in 0(\V |A 4 ) time by a call of Algorithm 1 with input G,V. 
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Algorithm 2 Compute vertex coordinates of H C U v< =wSv Q G 



INPUT: A graph G = (V,E), W C V; 

compute 0|s v (W)* and U ve wS v with Local Z>|s v (W')* computation and input G,W; 
H <- U V6W S V ; {Note W C V(H)}\ 
GoOn <— true 
while Go On do 

numjolass <— number of equivalence classes of 0| Sy (W)*\ 

Qi •(— subgraph of H induced by edges of <p ( - c X)\ Sy (W)* for all i = 1 to numjzlass; 
Qi(x) ■(— connected component of gi containing vertex x for each x € for all / = 1 to 

numjclass; 

9: if exist ; and j with n V(g_/(je))| > 1 for some x 6 then 

10: combine <p; and <p 7 -, i.e., compute U (pj in 0[5 V (W)*; 
it else 

12: GoOn false; 
13: end if 
14: end while 

15: v'o <— arbitrary vertex of W; 

16: label each vertexx in each Qi(\'o) uniquely with U(x) 6 {1, . . . , |<2i(vo)|}; 

17: set coordinates c r (vo) = for all r = 1 , . . . , numjclass 

18: for every vertex x G 2,(vo) and for all i = 1 to numjclass do 

19: set coordinates c r (x) = for all r = 1 , . . . , numjclass and r ^ ;; 

20: set coordinates ci(x) = U(x)\ 

21: end for 

22: c/ max max xe ^( H ) tifl ( v , x) ; 

23: Lj <- G V(H) I c///(vo,x) = i} for i = 1, . . . af max 

24: for i = 2 to L max do 

25: for all x G Lj that have not obtained coordinates yet do 

26: if for all // G N H (x) that already obtained coordinates holds (x,u) G <p r for some fixed r 
then 

27: set coordinate c r (x) = l r (x) {l r {x) is unique unused label}; 

28: set coordinates c, (x) = c,(«) for all ; = 1, . . . , numjclass, i ^ r; 

29: else if for all a G N H (x) holds u has not obtained coordinates then 
30: remove x and all edges adjacent to x from H; 

3 1 : remove x from L, ; 

32: else 

33: {now there are distinct neighbors u, w G N (x) and thus, have not been removed from H, 

such that they already obtained coordinates with ((x,u), (x,w)) g" 0|s r (W)*, i.e., (x,u) G 
9r, (x,w) G (p s , r^s} 

34: set coordinate c r (x) = c r (w); set coordinate c s (x) = Ci(w); 

35: set coordinates c,(x) = c,(m) for all ; = 1 to numjzlass, i ^ r,s; 

36: end if 

37: call ConsistencyCheck forx and vertices that already obtained coordinates; 
38: end for 
39: end for 

40: {H has been modified via deleting vertices x that fail the consistency checks.} 
41: OUTPUT: H with coordinatized vertices;; 
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Algorithm 3 ConsistencyCheck 

1: REQUIRE: Call ConsistencyCheck for vertex x from Algorithm 2; 

2: ENSURE: no two vertices obtain identical coordinates and adjacent vertices differ in exactly 
one coordinate; 

3: for all y e V(H), x^y that already obtained coordinates do 

4: {consistency check that no two vertices obtain the same coordinates} 

5: if c r (x) = c r (y) for all r — 1 to numjclass then 

6: remove x and all edges adjacent to x from H; 

7: remove x from L,; 

8: break for loop; 

9: else 

10: {consistency check that two adjacent vertices differ only in one r-th coordinate} 

11: if (x,y) is edge contained in some <p r and c r (x) = c r (y) or C;(x) ^ c,(y) for some i = 1 to 

numjclass, i ^ r then 

12: remove edge (x,y) from H; 

13: break for loop; 

14: end if 

15: end if 

16: end for 



As mentioned before, a vertex x of a Cartesian product Q" =1 G, is properly "coordinatized" 
by the vector c(x) := (ci (x), . . . ,c„(x)), whose entries are the vertices Cj(x) of its factor graphs G,-. 
Two adjacent vertices in a Cartesian product graph differ in exactly one coordinate. Furthermore, 
the coordinatization of a product is equivalent to an edge coloring of G in that edges (x,y) share the 
same color if x and y differ in the coordinate k. This colors the edges of G (with respect to the 
given product representation). 

Conversely, the idea of Algorithm 2 is to compute vertex coordinates of a subgraph of U v ewS v 
based on its 0| Sv (W)*-coloring. In particular, we want to compute coordinates that reflect parts of the 
5| Sv (W)*-coloring of U ve wS v in a consistent way. Consistent means that all adjacent vertices u and 
v with (u, v) G <p, c f|s,,(W0* differ exactly in their r-th position of their coordinate vectors, and no 
two distinct vertices obtain the same coordinate. This goal cannot always be achieved for all vertices 
contained in U vG iyS v . In [7, p. 280 et seqq.] a way is shown how to avoid those inconsistencies. In 
this approach colors of edges with "inconsistent" vertices are merged to one color. However, if the 
graph under investigation is only slightly perturbed, but prime, this approach would merge all colors 
to one. This is what we want to avoid. Instead of merging colors and hence, in order to preserve a 
possibly underlying product structure, we remove those vertices in U v6 iy5 v where consistency fails. 
This leads to a subgraph H C U v6 wS v where the edges are still D| Si (W)*-colored w.r.t. U v gjyS v and 
have the desired coordinates. In Algorithm 4 we finally compute Hi based on these coordinates and 
the edges of <p, c (D|s v (W)*)i#, 1 < i < k. Hence, the connected component of H induced by the 
edges of <p, e D\s v (W)* are subgraphs of layers Hj of the Cartesian product D^Lji/,- and therefore, H 
can be embedded into O k i=l Hi. 

Lemma 4.3. Given a graph G = (V,E) and ffCV such that (W) is connected, then Algorithm 2 
computes the coordinates of a subgraph H C G with H C U v( zwS v such that 

1. no two vertices ofH are assigned identical coordinates and 

2. adjacent vertices x andy with (x,y) S cp r - f \s y (W)* differ exactly in the r-th coordinate. 
The time complexity of Algorithm 2 is 0(\V | 2 A 2 ). 
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Proof. The init steps (Line 2 - 16) include the computation of i)|s v (W)*, H = U ve wS v , and the con- 
nected components Qi(x) that contain vertex x and which are induced by edges of (j?, e d\s v {W)*. 
By merging equivalence classes (Line 10) we ensure that after the first while-loop connected com- 
ponents induced by D|s v (W)* equivalence classes intersect in at most one vertex. Hence, vertices 
x in Qi(vo) can be assigned a unique label U(x) for each i = 1,. . .nwnxlass. In Line 17-21 we 
assign coordinates to each vertex contained in (2i( y o) for each i = 1, . . .nurruolass. Since any two 
distinct subgraphs 2;'( y o) and Qj(\'o) intersect only in vertex vo we can ensure that adjacent vertices 
in each subgraph <2i(vo) differ exactly in the z'-th position of their coordinate. We finally compute 
the distances from vo to all other vertices in H, and distance levels L, containing all vertices x with 
dn(vo,x) = i (Line 22 and 23). Notice, the preceding procedure assigns coordinates to all vertices of 
distance level L\. 

In Line 24 we scan all vertices in breadth-first search order w.r.t. to the root vq, beginning with 
vertices in Lq, and assign coordinates to them. This is iteratively done for all vertices in level L, 
which either obtain coordinates based on the coordinates of adjacent vertices or are removed from 
graph// and level L\. In particular, in the subroutine ConsistencyCheck (Algorithm 3) we might 
also delete vertices and therefore we have to consider three cases. 

First Case (Line 26): We assume that all neighbors of a chosen vertex x G L, that already obtained 
coordinates are contained in the same subgraph Q r (x). Hence, the coordinates of x should differ from 
their neighbor's coordinates in the r-th position. This is achieved by setting c r (x) to the unique label 
l r (x) and the rest of its coordinates identical to its neighbors. 

Second Case (Line 29): It might happen that vertex x does not have any neighbor with assigned 
coordinates, that is, either those neighbors of x are removed from H and Lj, j < i in some previous 
step, or they have not obtained coordinates so far. If this case occurs, then we also remove vertex x 
from H and L,, since no information to coordinatize vertex x can be inferred from its neighbors. 
Third Case (Line 32): Let u,w € N (x) be neighbors of x such that u and w have already assigned 
coordinates and the edges (x,u) and (x,w) are in different equivalence classes. Assume (x,u) G (p r 
and (x, w) G <p s , r/i. Keep in mind that x should then differ from u and w in the r-th and in the s-th 
position of its coordinates, respectively. Thus, we set coordinate c r (x) = c r (w) and c s (x) = c s (u). 
The remaining coordinates of x are chosen to be identical to the coordinates of u. Note, we basically 
follow in this case the strategy to coordinatize vertices as proposed in [1]. 

In order to ensure that no two vertices obtained the same coordinates or that two adjacent 
vertices differ in exactly one coordinate we provide a consistency check in Line 37 and Algorithm 
3. If x has the same coordinate as some previous coordinatized vertex we remove x from H and L,. 
If x has a neighbor y with coordinates that differ in more than one position from the coordinates of x 
we delete the edge (x,y) from H. 

To summarize, we end up with a subgraph H C \J ve wS v , such that the vertices of H are uniquely 
coordinatized and such that adjacent edges (x,y) G q> r - Q\s r (W)* differ exactly in the r-th position 
of their coordinates. 

We complete the proof by determining the time complexity of Algorithm 2. Lemma 4. 1 implies 
that Algorithm 1 determines 0i Sv (W)* and U ve wS v in 0( \V |A 4 ) time. Since (W) is connected, Lemma 
3.7 implies that D|s v (W)* has at most A equivalence classes and therefore, the while-loop (Line 5 - 
14) runs at most A times. The computation of the graphs Q, and Qi(x) within this while-loop can be 
done via a breadth-first search in 0(\E\ + \V\) = 0(|V|A) time, since there are at most | V [A edges 
and connected components. The intersection and the union of Qj and Qj can be computed in 0(\V | 2 ). 
Hence, the overall-time complexity of the while-loop is <9(A|V| 2 ). The assignments of coordinates 
to vertices x G 2;(vo) can be done in 0(A) time. Since there are at most |V| vertices and at most 
A equivalence classes we end in (9(|V|A 2 ) time. Computing distances from vo to all other vertices 
and the computation of L,- can be achieved via breadth-first search in 0(\E\ + \V |) = 0(|V|A) time. 
Consider now the two for-loops in Line 24 and 25. Each vertex is traversed exactly once. Hence 
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Algorithm 4 Embedding of H into Cartesian product 

1: INPUT: A graph G = (V,E) with coordinatized vertices; 

2: for each position i = 1 to r of coordinates do 

3: initialize graph //, = 0; 

4: for each vertex v G V do 

5: if c/(v) {Hi) then 

6: add Cj(v) toV(#,-); 

7: end if 

8: end for 

9: end for 

10: for each position ; = 1 to r of coordinates do 

11: for each edge (x,y) G E do 

12: if c;(x) 7^ Cj(y) and edge (c,-(ji;),Cj(y)) ^ £(//,■) then 

13: add (c/(ac),Ci(y)) to £(#,-); 

14: end if 

15: end for 

16: end for 

17: OUTPUT: Factors //, and Cartesian product U r j=l Hi where G can be embedded into; 



these for-loops run 0(|V|) times. For each vertex in each distance levels we check whether there 
are neighbors in level Lj_i, which are at most A for each vertex x, and compute the A positions of 
the coordinates for each such vertex. The consistency check (Algorithm 3) runs in 0(|V|(A + A) = 
0(|V|A) time. Hence, the overall time complexity of the for-loop (Line 24 - Line 39) is 0(|V| 2 A 2 ). 

Combining these results, one can conclude that the time complexity of Algorithm 2 is 
0(|V| 2 A 2 ). □ 

Lemma 4.4. Given a graph G = (V,E) with bounded maximum degree A obtained from Algorithm 
2 with coordinatized vertices. Then Algorithm 4 computes factors Hi such that G can be embedded 
into W i=l Hi in 0(|£|A) time. 

Proof. After running Algorithm 2 we obtain a graph G = (V,E) such that vertices x G V have con- 
sistent coordinates c(x) = (c\ (x), . . . ,c r (x)), i.e, no two vertices of G have identical coordinates and 
adjacent vertices x and y with (x,y) G (f>, differ only in the ;-th position of their coordinates. 

We first compute empty graphs H\,...,H, and add for each vertex x and for each Cj(x) of 
its coordinates c(x) = (c\(x), . . .c,-(x)) the vertex e;(x) to Hi. Different vertices c,(x) and Ci(y) are 
connected in Hj whenever there is an edge (x,y) G E. We define a map y : V(G) — > V(H) with 
x M> c(x). Since no two vertices of G have identical coordinates y is injective. Furthermore, since 
adjacent vertices x and y that differ only in one, say the ;-th, position of their coordinates are mapped 
to the edge (c ; (x),c ( (y)) contained in factor //, and by definition of the Cartesian product we can 
conclude that the map y is an homomorphism and hence, an embedding of G into H. 

The first two for-loops run \V |A times, that is 0(|,E|). The second two for-loops run |£|A times, 
hence we end in overall time complexity of 0(|£'|A). □ 

To complete the paper, we explain how the last algorithms, in particular, Algorithm 1 , 2 and 
4 can be used as suitable heuristics to find approximate products; see also Figures 5 and 6. Note, 
by Corollary 4.2 Algorithm 1 can be used to compute 8(G)*. However, most graphs are prime and 
8(G)* would consist only of one equivalence class. Thus we are interested in subsets of 8(G)* which 
provide enough information of large factorizable or "into non-trivial Cartesian product embeddable" 
subgraphs. This can be achieved by ignoring regions S v where 5| 5v has only one or less than a 
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(a) A Cartesian prime graph G = (V,E) is shown. For all 
vertices x EV (marked with "X") the respective D\ Ss has 
only one equivalence class. Thus, we use only all non- 
"X"-marked vertices, pooled in the set W C V and call 
Local 0|5',.(W)* computation (Alg. 1). The equiva- 
lence classes of 0| Sr for vertex v = vo are highlighted by 
dashed and thick edges. 
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(c) Shown is the graph G with coordinatized vertices for 
all x £ UjLjL,. Note, the vertex x with coordinates (37) ob- 
tained a new unused second coordinate 7, since all edges 
(u,x) where u already obtained coordiantes are from the 
same equivalence class (Alg. 2, Line 26). Thus, coordi- 
nates cannot be combined. 
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(b) After calling Local 0|s r (W)* computation 
(Alg. 1) we obtain the equivalence classes of 0|s v (W)* 
highlighted by dashed and thick edges. After calling 
Compute vertex coordinates (Alg. 2, Line 15 - 
21) we obtain a graph where the vertices in each GJ -layer 
obtain unique coordinates. 
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(d) Shown is the graph G with coordinatized 
vertices for all x 6 U^ =1 L,. Note, after running 
ConsistencyCheck (Alg. 3, Line 11) the edge 
between the vertices with coordinates (37) and (25) 
is deleted, since the vertices differ in more than one 
coordinate. 



Figure 5. The basic steps of Algorithm 1 and 2 



given threshold number of equivalence classes. Hence, only subsets W CV where 0|s v (W)* has a 
sufficiently large number of equivalence classes are of interest. For this, we would cover a graph 
by starting at some vertex v <E V, compute S v and d\ Sv , and check if J)i Sv has the desired number of 
equivalence classes; see Figure 5(a). If not, we take another vertex w € V and repeat this procedure 
with w. If D|5 v has the desired number of equivalence classes we would take a neighbor w of v, 
compute S w and 0| Su , and check whether (0| Sh , U Q\s y )* has the desired number of equivalence classes. 
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Figure 6. After running Algorithm 1 and 2 we obtain H as a subgraph of the 
graph G in Figure 5, with coordinatized vertices, and edges colored w.r.t. f|s r (W)* 
equivalence classes. After running Embedding of H into Cartesian 
product (Alg. 4) we obtain the putative factors Hi and H2 of H and, hence, 
of G. Note, due to the coordinatization of H the embedding of H into H\ UH2 can 
easily be determined. 

If so, then we continue with neighbors of v and w and to extend the regions that can be embedded 
into a Cartesian products. To find such regions one can easily adapt Algorithms 2 and 4. 

Note, after running Algorithm 1 one could take out one of largest connected component of 
each equivalence class induced by edges with the respective "colors" to obtain putative factors; see 
Figure 5(b). However, even knowing putative factors does not yield information about which edges 
should be added or deleted to obtain a product graph. For this, coordinates are necessary. They can 
be computed by Algorithm 2 and used as input for Algorithm 4; see Figure 6. 

Finally, even the most general methods for computing approximate strong products only com- 
pute a (partial) product coloring of the graphs G under investigation. They yield putative factors, 
but no coordinatization [8]. However, Algorithm 4 can be adapted to find the coordinates of the so- 
called underlying approximate Cartesian skeleton of such graphs, and can thus be used to find an 
embedding of (the approximate strong product) G into a non-trivial strong product graph. 
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